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Abstract. We study stimulated light scattering off a superfluid Fermi gas of atoms 
at finite temperature. We derive response function that takes into account vertex 
correction due to final state interactions; and analyze finite temperature effects on 
collective and quasiparticle excitations of a uniform superfiuid Fermi gas. Light 
polarization is shown to play an important role in excitations. Our results suggest 
that it is possible to excite Bogoliubov- Anderson phonon at a large scattering length 
by light scattering. 
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1. Introduction 

Since the first realization of Bose-Einstein condensation in atomic gases in 1995 pLj, tliere 
lias been tremendous growth in research activities with cold atoms. Recent experimental 
studies [21 IHl m [HI m IZl [H] with cold fermionic atoms have generated renewed interest in 
quantum many-body physics. Atomic Fermi gases in traps provide a unique laboratory 
system for exploring physics of interacting fermions with tunable interactions. Fermi 
degeneracy in a trapped atomic gas was first demonstrated by DeMacro and Jin in 
1999 p. In recent past, there has been many reports of possible observations of Fermi 
superfluidity (FS). However, an unambiguous evidence of FS was shown by Zwierlein et 
al. yiOj. The detection of pairing gap [mini [131 [H [15] and collective modes [I6l[171[l8] 
of FS are current issues of interest. Physics at the crossover [191 [201 [21] between BCS 
state of atoms and EEC of molecules [22] formed from Fermi atoms is of prime interest. 
A number of theoretical investigations [23] have dealt with FS near crossover. Two 
very recent experiments [211 EH] on two-component Fermi gases with imbalanced spin 
components provide new insight on the nature of FS and perhaps indicate to some new 
state of matter known as interior gap (IG) superfluidity predicted by Liu and Wilczek 
[26j . The occurrence of IG superfluidity in a two-component Fermi gas was theoretically 
predicted by Deb et al. [27]. 

In order to study the nature of FS, it is important to derive the appropriate response 
function of Cooper-paired Fermi atoms due to an external perturbation (such as photon 
or rf field). Our purpose here is to calculate response function of superfluid Fermi gas at 
finite temperature due to stimulated light scattering. In a previous paper [28] , we derived 
response function at zero temperature. We have also shown that it is possible to excite 
selectively a single partner atom (of a particular hyperfine spin state) of a Cooper-pair 
exploiting light polarizations in the presence of a strong magnetic field. We present here 
detailed method of calculation of response function at finite temperature due to light 
scattering. We study the effects of finite temperature and light polarization on single- 
particle excitation as well as collective mode of density fluctuations. This collective mode 
known as Bogoliubov- Anderson (BA) phonon [29j has been theoretically studied in the 
context of fermionic atoms [30] . We here show that it may be possible to excite this mode 
by light scattering: This mode appears as a resonance in dynamic structure versus energy 
transfer in scattering. At finite temperature the resonance becomes broadened due to 
Landau damping. In case of single-particle excitation spectrum when Cooper-pairs are 
broken {uo > 2A), the peak occurs at a higher energy as the temperature is lowered. 
If the momentum transfer is higher, single-particle spectrum becomes sharper. Light 
polarization significantly affects the excitation spectrum. Circular polarization of light 
leads to a positive shift of the peak of single-particle spectrum at lower temperatures. 

The paper is organized in the following way. In the following section, we discuss in 
brief the stimulated light scattering off Cooper-paired Fermi atoms. In Sec. II, we present 
in detail the derivation of response function with vertex correction at finite temperature. 
We next describe numerical results in Sec. III. The paper is concluded in Sec. IV. 



Finite temperature effects in light scattering off Cooper-paired Fermi atoms 



3 



2. Stimulated photon scattering 

Stimulated light scattering will occur when two laser beams with nearly equal frequencies 
are impinged on atomic gas and tuned far-off the resonance of a transition frequency of 
the atoms. Let us specifically consider two-component Fermi gas of ^Li. Here the two 
components imply the two hyperfine spin states of F = 1/2 ground level. As discussed 
previously [28] , large Zeeman shifts of hyperfine sub-levels near Feshbach resonance allow 
us to utilize light polarization to control scattering to a significant extent. Fig.l of Ref. 
[28] describes polarization-selective light scattering by which the amount of momentum 
and energy transfer to the two partner atoms of a Cooper-pair can be controlled. It 
is thereby possible to scatter photons from either spin component, however because 
of long-range correlation between two components due to s— wave Cooper-pairing, the 
other spin component will also be affected by photon scattering. When light fields are 
treated classically, the effective atom-field interaction Hamiltonian is 

where ao- k(oCTk) represents the annihilation(creation) operator of an atom with hyperfine 
spin cr and center-of-mass momentum k, q is the momentum transfer and 70-0- is the bare 
vertex corresponding scattering without any change in spin state a. 

To describe density-density correlation function, we define the density operators by 

= Ea,k al,k+q«a,k and 

Pq^ = 7<x<7al,k+qa^,k (2) 

k,rj 

3. Response function 

The response function due to density fiuctuation is 

x(q,r-r') = -(T.[p(^)(r)p^"i(r')]), (3) 

where (■ ■ ■) implies thermal averaging and Tr is the complex time r ordering operator. 
By generalized fiuctuation-dissipation theorem, the dynamic structure factor is given by 

S{q,uj) = --[1 + nB{uj)]lm[x{q,z = u + 16)], (4) 
vr 

where x(ct) -2) represents Fourier transform of 
3.1. Vertex equation 

Within the framework of Nambu-Gorkov formalism of superconductivity, using Pauli 
matrices Tj, the vertex function can be written as [31] 

r(fc+,fc_) = 7- j r'^E^3G(A:;)r(fc;,fc^)G(fc'_)r3r(k,k'), (5) 

where k+ = {k + q/2,i(ci;„ -|- z/m/2)} and k_ = {k — q/2,i{ujn — ^'m/2)} with 
^n/fi = (2^ + and z^m/^ = '2m{hj3)~^ representing the fermionic and bosonic 
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Matsubara frequency, respectively. The Green function can be expressed in a matrix 
form as 

"^'^ = A^l • 



where /i„/^ is Matsubara frequency, = + A| and = — 1^- with = 
fi^ k'^ / {2m) . The bare vertex can be written in a matrix form 

7 = 70^0 + 73i"3, (7) 

where 70 = [7tt - 7ii]/2 and 73 = [7tt + 7ii]/2. 

The response function is related to the vertex function by 

X(q,^) = / ^r^ETr[7G(k+)r(k+,k_)G(k_)] (8) 

To solve the vertex equation, let us expand the vertex function in terms of Pauli 
matrices as 

r(A;+,fc_)=f]r«(k,q,iz/Jr,. (9) 

We replace l^(k, k') by an effective mean field potential K//. In the weak-coupling 
limit, this potential reduces to the form Kveak = QO-s where g = Anfi^ /m and is 
s-wave scattering length. Since we are interested only in Cooper-pairing regime, we 
consider attractive interaction only and hence is assumed to be negative. On replacing 
l^(k, k') by gosi the gap A becomes fc- independent. The vertex function also becomes 
fc— independent but remains a function of q and ivm only. 
Using Eqs. ([6]) and ([9]) in Eq. ([5]), we can write 

r(k, q, lUm) =lk- gas J 

-1 1 ^ 

where ^k = h^k'^ / (2m) — yU and 

+ + Z/^/2) + l{uJn - yj2)]V^^^ 

+ [ik'^ikL - K + vJ1){yJn - vjl) + A2]r(°), (11) 

1/1 = - %\%{uj^ + z/^/2)efc^ - z(^„ - z/„/2)efc;]r(2) - A[efc; + ^fc'jr^^) 

- A[z(a;„ + Vml2) + 2(^„ - z/^/2)]r(o) 

+ \ik'^ik'_ + (^n + ymm{^n ' ^^m/2) - A^]!^!), (12) 

2/2 = - ^A[z(o;„ + z/„/2) - «(cu„ - z/™/2)]r(=^) - i^\ik'^ - e^jr^") 

- z[4/^i(u;„ - z/^/2) - ii(u;„ + z/„/2)^fc--']r^^^ 

+ [e^Vefc'. + (^n + Vrnmi^n ' ^^m/2) + A^]!^^), (13) 
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There are basically two types of Matsubara frequency sums: 



7i(k',q) = /3-i^ 

n 

/,(k',q)=r^^ 



1 



{uJn - 1^J2)^ + El_ 
1 



X 



{un + Vml2Y + - i/^/2)2 + £;2^ 

Now, the term like uj/[uj'^ + can be written in a separable form 



1 

2i 



+ - 



iuj + E^. iuj — El 



while the term hke + El] can be decomposed as 



u^ + El 2Ek 



(14) 

(15) 
(16) 

(17) 
(18) 



.iu + Ek icu — Ek. 

The terms which are odd in frequency will not contribute to the sum and so can be 
omitted. Using these decompositions and after some algebra as done in Appendix- A, 
we can express 



/i(k',q) 
/2(k',q) 



-[Tn + ri2 + r2i + 



1 



11 



■ 12 



^21 



where, 



Til 
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^21 = 



Too — 



4Ek'^Ek'_ 

tanU 3 Ei,'_ /-I) - taiili(;7£/,;/2) 
2(£;fc; - Ek'_ + ium) 

tanh(/3Efc/_/2) + tanh(/^Efc^/2) 

2{Ek'^ + Ek'_ + tUm) 

tanh(/3£;fe^/2) +tanh(/3£;fe/_/2) 
2{Ek'^ + Ek'_ - ivm) 

tanh(/5Efc/_/2) - tanh(/?Efe/ /2) 



^22 



(19) 
(20) 

(21) 
(22) 
(23) 
(24) 



2(-Efc^ - Ey_ - ip„ 

In what follows, we use as the unit of energy the Fermi energy eF = h'^kp/ (2m) and 
accordingly scale all the quantities. We denote A = A/ep, = ^k/^F, k = k/kp and so 
on. Let X = qk cos 9^ where 9 is the angle between k and q. For notational convenience, 
let £ = lk + gV4, ^1 



Ek 



{S - xy + A2 and E2 = E, 



k+ 



{s + xy + A^. 
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Having now performed the Matsubara frequency sum, omitting the terms which are 
odd in X, we can express the vertex equation as 

r d^k' ^ 

r(k, q, ium) =1 / TTT-T^ Ci(k')ri, (25) 



where = goskp/ep and 

02 = Az/„/2r(3) - u^gi^r^'-) + [{s' -x' + A')h + /i]r(^), 
C3 = 2A^/2r(i) + Az/„/2r(2) + [{s^ -x^- A^)h - /i]r(=^) 

Now, the vertex terms T^^^ form four coupled algebraic equations 

r d^k' 

T^'^ =7^-'^sJ J^C,, 2 = 0,3 (26) 
r d^k' 

In the limit g ^ 0, and T^"^^ A, the equation for T^"^^ reduces to the standard 

BCS gap equation 

~ r d^k tanh((3EJ2) ~ , , 

A = -Ks 7 — ^ ' ' a 28 

i (27r)3 2Efc ^ ' 

which implies |ks|/o = 1 where 

. _ [ d^k tanh(/3Efc/2) 

The gap equation expressed in this form has logarithmic divergence. However, this 
divergence can be removed by renormalizing the mean-field interaction via subtracting 
the zero temperature and zero pairing-field (A = 0) part of the integral. Although this 
gap equation resembles the standard weak-couphng BCS gap equation, the chemical 
potential /i can deviate significantly from its weak-coupling value /i ~ ei?. The chemical 
potential is given by single-spin superfluid number density 

n = \l f 1 - 1^ tanh(/5Efc/2)) . (30) 



2 y (27r)3 V J 

Let us first consider the off-diagonal terms F^^^ and F^^^. Making use of the relation 
f l29|) and the analytic continuation ium + iO^ , 

^ JiF(^) + 2A JiF(3) 

(J2-A2j2 + /l)-/o 

^(2) ^ 2A(^ J2F(3) - JiF(^) 

(/2 + AV2 + /l) - Jo 
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hiq^ii^m 
Ji(q, zz/„ 
J2{q.,iu„ 



(27r)3 
_fk_ 

d^k 



(27r)3 
(2^ 



/i(k,q) 
^/2(k,q) 
/2(k,q) 



(33) 
(34) 
(35) 
(36) 



The method of calculation of various integrals is described in Appendix-A. Eliminating 
from Eqs. (13T1) and (1321) . we have 



r(2) = tACu 



-1 



A 



2Jl 



r(3) 



(37) 



where D± = (Ji + I2 — Iq) ± A^J2. Here we note that although the integral /i, I2 and 
Jo have logarithmic divergence, this divergence does not pose any problem. Because, at 
the end all the divergences are exactly canceled out and thus all the vertex terms F^^^'s 
and the response function remain finite. The integrals Ji and J2 are finite. From Eq. 
fl26l). we have the solution 



r(o) 



7o 



1 + KsB 



where 



B = h + AV2 - h. 



On substitution of Eqs. (!3T|) and (!32|) . we have the solution 

r(3) ^3 



l + KsF 



(38) 
(39) 
(40) 



where 



A + 



X 



D,D_ 



n -1 



2{ACoJi] 

(2AJi 
D_ 



A 



_2Jl_ 
D.D. 



+ 



with 



(41) 



(42) 



A = l2- AV2 - h 

After having done Matsubara frequency sum and analytic continuation, we can 
write the response function as [see Appendix-B] 

X{q, 00) = 257oF(°) + 2A-^3T^^^ - 2iuKj2i^V^'^'^ + 4A Ji73F(i) (43) 
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Figure 1. Gap A (in unit of tp) is plotted as a function of temperature ksT (in unit 
of ep) for feFl^sl = 2 (a) and fcir-lasl = 4 (c), while subplots (b) and (d) exhibit the 
variation of the chemical potential /i as a function of temperature for fc_F|as| = 2 and 
fcj-jasl = 4, respectively. 



On replacing F*^^^ and T^^^ in terms of V'^'^\ we have 
x(g,u;) = 257orW + 2F73r(3) 



7o + (44) 



l + KsB"' l + KsF 

Now, the dynamic structure factor can easily be written as 
5'(q,t^)= - -[1 + ni3(u;)]Im[x(q,u;)]. 

TT 



2 

-[1 +nB(tu)] 

TT 



7o'lm(B) ^ 7|lm(F) 



|1 + K,5|2 |1 + /€,F|2 



(45) 



4. results and discussions 



Before we elaborate our results, we note a few pertinent points. First, in calculating 
response function, unlike weak-coupling case, all the energy integrations are carried out 
over the entire energy range meaning that ^ ranges from — yU to oo. Second, as regards 
vertex correction in light scattering, our formalism of calculating response function 
can account for any arbitrary polarization of incident light. Third, we have devised 
a procedure whereby we carry out first the angular integrations by parts considering 
C(j as a complex parameter z as described in Appendix-A. This leads to trigonometric 
functions of which are then easily analytically continued using the limit z ^ uj + iO"'". 
This is unlike the usual approach where the functions like l/[u; ± {E±) + zO^] are first 
separated into real and imaginary parts. The imaginary part is a delta function, the 
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Figure 2. Subplot-(a): Dynamic structure factor S{Lo,q) (in dimensionless unit) is 
plotted as a function of energy transfer ld (in unit of tp) for q = Q.2)kp^ kp\as\ — 2.0, 
ksT — 0.009ei? (Tc = 0.0265ei?) for unpolarized light with 711 = 722 = 7, that is, 
7o and 73 — 7. Subplot-b: Same as in (a), but for circularly polarized light 
with 711 = and 722 7^ 0, that is, I70I = I73I = 722/2. Subplots (c) and (d) are the 
counterparts of (a) and (b), respectively; but for fc^T — 0.021. 



energy and angular integrations over which require finding the roots of a comphcated 
polynomial equation. 

For all our numerical illustrations, the coupling strength Kg is taken to be negative 
since we are interested in Cooper-pairing regime only. Figure 1 shows temperature 
variation of A and /i for two different scattering lengths a^. Comparing Fig.l (a) 
with Fig.l (c), we note that critical temperature becomes large for large scattering 
length. This implies that finite temperature effects are particularly important for large 
scattering lengths. In the zero temperature limit, the gap becomes almost independent 
of temperature, but remain largely dependent on interaction parameter. As in our 
previous paper [28j for T = 0, for numerical illustration we consider two cases: Case-I: 
When both the laser beams are unpolarized, that is, 7ii(= 7tt) = 722 (= 7u)j Case-II: 
When both the beams are either cr+ or cr_ polarized. These are two limiting cases only, 
we reemphasize that our formalism can handle any arbitrary polarization. 

Let us next discuss the regime of collective excitation of density-density fluctuations 
for energies uj < 2 A. This regime is characterized by long-wave mode of vibration known 
as Bogoliubov- Anderson (BA) phonon which results from vertex correction. The general 
expression for x as given by Eq. (H5ll has two parts which are proportional to 7o and 7I, 
respectively. Since the real part of B is always negative, the flrst part is always flnite. 
When T^^^ 70, F*^^) 73 and F* with i = 1,2, that is, when vertex correction is 
neglected, the response function reduces to 

X(g,u;) = 2572 + 2^72 (46) 
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Figure 3. S{uj, q) as a function of w for two different values of q as shown in the figure. 
The other parameters are kplusl = 4.0, ksT = 0.055eF (T^ = 0.1246;^, A = 0.2036^), 

711 = 722 = 1- 




Figure 4. Quasiparticle excitation spectrum (w > 2A) for a uniform Fermi superfluid 
at two different temperatures fc^T = 0.115eF (a, b) and ksT = 0.075eF (c, d), but 
for the same A;F|as| = 4.0. The subplots (a) and (c) are for unpolarized Hght with 
7ii = 722 = 7, while (b) and (d) are for circularly polarized light with 711 — 0. The 
dashed curve refers to q = OAkp and the solid one to q = Q.Gkp in all the subplots. A 
is 0.098ei? and O.ISScf at temperatures O.llbep and 0.075eF, respectively. 
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which has no pole. However, when vertex correction is added, there arises the pole of x 
which is given by 

1 



(47) 

I '^s I 

In the limit \as\ — > 0, the pole is determined by = 0. The real part of the pole will 
correspond to the phonon energy while the negative imaginary part will give the damping 
of the mode. Only in the low g-regime and T < T^, BA mode will be well defined. For 
larger temperatures (but < Tc), this mode will be ill defined because of large Landau 
damping which occurs due to its coupling with thermally excited quasiparticles. At zero 
temperature, in the limit — ^ 0^, we have reproduced the standard result 

t^BA = -^hqvF (48) 

as calculated in Appendix-C. The BA phonon appears as a resonance in the spectrum 
of dynamic structure factor as is demonstrated in Figs. 2 and 3. Comparing Fig. 2(a) 
with Fig. 2(c), we notice that as the temperature is increased, the spectrum becomes 
broadened. The resonance energy does not depend on the state of polarization of light. 
However, the peak and width of the resonance may depend on light polarization. Figure 
4 displays the quasiparticle excitation spectrum when uj exceeds the pair-breaking energy 
2 A. As temperature is decreased, the peak makes a positive shift because of increase in 
A. For larger q values, the peak is larger. At lower temperatures, polarization of light 
significantly affect the spectrum as can be noticed by comparing Fig. 4(c) with 4(d). 



5. conclusion 



In conclusion, we have derived finite temperature response function of superfluid Fermi 
gas due to scattering of polarized light. The response function takes into account vertex 
correction due to final state interactions. We have presented selective results on dynamic 
structure factor (DSF) deduced from the response function with the aid of generalized 
fluctuation-dissipation theorem. In-gap collective mode of density fluctuation, known 
as Bogoliubov-Anderson phonon, appears as a strong resonance in the plot of DSF 
as a function of energy transfer uj at low momentum transfer q. As the temperature 
increases, the width of the resonance increases due to Landau damping. At a large 
scattering length (/ciT'losI > 1), the resonance may occur at a finite and appreciable 
value of q. Polarization-selective light scattering may be useful in exciting BA mode. 
We have also presented results on single-particle excitation spectrum when Cooper- 
pairs are broken {uj > 2 A). As the temperature is decreased, the peak of the spectrum 
makes a positive shift due to increase of pair-breaking energy (2A). We have also 
shown the effect of light polarization on collective as well as single-particle excitations. 
Particularly important is the positive shift of the peak of the single-particle excitation 
spectrum at lower temperatures due to circularly polarized light as compered to that 
due to unpolarized light (Fig. 4). Furthermore, higher the momentum transfer, sharper 
is the single-particle excitation spectrum. As pointed out earher [27], in stimulated 
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light scattering, the momentum transfer may be as large as comparable with the Fermi 
momentum. Since in polarization-selective stimulated light scattering, only the atoms 
with the same spin state are predominantly scattered [28], this kind of light scattering 
with large momentum transfer may be useful in detecting the gap energy and its 
temperature dependence. In the present paper, our study is confined to uniform FS 
only. Local density approximation (LDA) may be applied for studying single-particle 
excitation of a trapped FS at a large momentum transfer as done in [28]. However, for 
studying collective mode of trapped system, LDA may be a bad approximation. We 
hope to address this problem in our future communication. 
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Appendix-A 



Here we first calculate the four Tij terms as expressed in Eqs. ( l2niMl) and then describe 
the method of calculation of various integrals. We can express 

1 1 



Similarly, 



11 



12 



'21 



T, 



— y 

Ek'^ - Ey_ +iymn 

Hpi-Ek'J - npi-Ek'^) 
Ek'^ — Ek'_ + ium 

npiEk'J - npi-Ek'^) 
Ek'^ + Eki_ + ivm 

npi-Ek' ) - npiEk'J 

—Ey^ — Ek'_ + ih'm 
npiEk'^) - upi-Ek'J 
Ek'^ + Ek'_ — ium 

npiEk') - upiEy ) 



1 



k' 



i{iUn - + Ek'^ i{uJn + iym/2) + Ek' 



(A.l) 



(A.2) 



22 



Ey — Ey 



(A.3) 



(A.4) 



One can write np^E) = l/2-tanh(/3E/2)/2 and npi-E) = l/2 + tanh(/3E/2)/2. Using 
these relations, we obtain the Eqs. fl2m24p . 

Let us now consider the integral f and Jj. Setting z = ivmi we can write 

1 



i2ny 



sin 6d9k dkli{k, q, z) 
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where 



Here 



J2{q, z) = ^ / edk £ dy^n. {k, q, y, z) 
n^k, q, y, z) = + ± ^) - 



T± = tanh(/5E2/2) ± tanh(/3Ei2), 
Here y — cos ^. Similarly, we can express 



-^2(9,^) 



4(27r)^ 



1 - x"^ 
k'^dk I dy 7^_ (/c, g, |/, 2;) 

— 1 il/2-C/i 



14 
(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 
(A.IO) 



In the limit T ^ 0, we have = 2 and T_ = 0. Now, 



k dk I dy 



X 



4(27r)2 

'T+{Ei + E2){{E2 - ^i)2 - + T_{E2 - E,){{Ei + E^)^ - z^} 
{El - Ely - 2z\El + El) + z^ 
1 r 2 

J ^'"^^ 7-1 ^y{El-Efy-2z^{E! ' ^2 



X 



£^2 tanh 



Substituting E^ — Ef — A£x and Ef + E^ — 2£^ + 2x^ + 2A'^, the above equation can 
then be expressed as 



3po 



J J—gk 



dx- 



fk z^{z^ - - 4A2) + 4(4£:2 - z^)x^ 



X 



-z^J^. - ^ExT- 



f 1 f'^i 1 



where 



3po 

^2 - 4g2 
z\S^ - ZD 



(A.ll) 

(A.12) 
(A.13) 



= z2j^^+ + ^ExT-, (A. 14) 

with jFj_ = Ex tanh(/?_Ei/2)±i?2 tanh(/3i?2/2). Here po = ^f/ (67r2) is the number density 
of noninteracting Fermi gas of a single-spin. We further note that E\{x) = E2{—x) 
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implying J-±{x) = ±J-'±{—x). Let us carry out x— integration by parts, assuming 
and (1 — u'^x'^)~^ be the first and second integrand, respectively. Then we have 



rqk 

Jo 



dx ^(x 



1 



X 



1 + i/^X^ 

rqk 

^( qk) tein~^ ( uqk) - / dx^' (x) tan"^ (ux) 
Jo 



(A.15) 

Let us now consider the analytic continuation by taking the limit z —> uj + iO^, 
where Co is the energy transfer. There are two regimes of excitations: I. Quasi-particle 
regime: a) > 2A and IL Collective oscillation regime: oj < 2A. We first concentrate 
on the regime of quasi-particle excitation, that is, Cu > 2A. In this regime. So > 0. 
Let (j) — tan~^(z/a;) and So — \Jcj'^ — aA?/2. We have three cases for consideration of 
analytic continuation: (1) For \S\ < So, we have 

1 



IVi — i- 



4£:2 



itanh ^{uix), u^x < 1, 

— hitanh"^ ( ^ , Pix > 1, 

2 \vixj 



(2) For So < \S\ < Cj/2, we have 



1 cD2-4£:2 

tan0 — > V2X, 
(3) For \S\ > Cj/2, we have 



1 /4<f2 



~2 



(j) — > itanh~"^(i/3a;), v^x < 1, 



— > — hi tanh . u^x > 1, 

2 V 1^3.7;/ 



(A. 16) 
(A.17) 

(A. 18) 

(A.19) 
(A.20) 

(A.21) 
(A.22) 

(A.23) 



After having done angular integration and analytic continuation, the energy 
integration is carried out numerically. Let us next consider the integral J2. Towards 
this end, we have 



7^_ 



z'^J+{x) +4:SxJ4x) 



E1E2 z^{S^ -Si){l + u^x^) 
where J7±(a;) = tanh[/5i?2/2]/£'2 ± tanh[/5£'i/2]/£^i. Using this, we can write 

where 



' IQtpqJ z\S^-Z^)Jo 
^!{x) = z'^J+ + 4SxJ_ 



(A.24) 

(A.25) 
(A.26) 
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Thus J2 integral can be evaluated exactly the same way as done for Ii. Similarly, we 
can find the other two integrals I2 and Ji. 

Appendix-B 

We can rewrite the response equation 

^^^'^^^ /(2;0^^ '?[(a;„ + ^.„)2 + £;|]K + £;?] 

X IV [7/C{i(u;n - z^m/2)ro + ^k.ra + An}] (B.l) 

where 

j=0 [ j=l,2 

+ A (r(°)ri + r(^)ro + ^ E r(^')ei,r„^. 1 (B.2) 
where rij = |1 + €ijj\ and €ij — —eji = 1 if is (1,2) or (2,3) or (3,1). Now, 

K, X - i^m/2)To + Ck-Ts + Ati} = i{Un - Vm/'^)^ 

+ ikAi^n + i^m/2)) I r(o)T3 + r(3)ro + ^ ^ e.ar^^Va-, 

[ i=i,2 

+ 6,6_|r(°Vo + r(^V3- E r^^'^^.j 

+ r(^)T3 - zr(2)ro + r(3)ri} + jr^^Vo + r^^Vi - E r^^V.j (b.3) 

Tr [7/C(i(a;n - i^m/2)ro + ^fe_r3 + An)] 
= -2(a;„-i/„/2)(^„ + i/^/2) ^ 7,r(^'^ 

j=0,3 

+ 2^(o;„ - z/^/2)a+[7or(=^) + 73r(°)] + 2i(u;„ - v^/2)K[^oT^^^ + i73r(2)] 
+ 2efe_^((a;„ + P^I2)) {731(0) ^ ^^r(3)| 

+ 2^+^- {r(°)7o + r(-^)73} + 2Aa_ {r(i)73 + i^^^ho] 
+ 2iA(^„ + vj2) {r(i)7o - iv^^h^} 

+ 2Ae,^ {r(^)73 - ir(^)7o} + 2A^ {r(o)7o - v^'^n} (b.4) 



So, 
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Appendix- C 



Here we turn our attention to collective oscillation regime where Sq < 0. We do not make 
any attempt to evaluate analytically the collective mode for any interaction strength 
(/ci?|as|) at finite temperature. However, it is possible to calculate analytically BA mode 
energy in the weak-coupling limit /cfIosI — > at zero temperature. Towards this end, 
let us first calculate D^{z — > cu + iO'^). Restricting a; < 1, we expand the functions 
E± = £'2,1 (x), and 1/E±{x) up to second order in x 

± 2£x £^x^ 







1 


1 




~ E 



2E 



2E^ 



x^±2Ex 3£:V 

1 r^3^ + 



At zero temperature, we obtain 



2E^ 



$(x ~ 0) ~ i:?{2E + — ) - 8£:^^ - 

E E 



x^ uP'E'^x^ 



x^ 



x 



E3 ' 



We can write I1 + I2 + J2 

3po 



4e, 



X 



E^ 

1^ + I^^, where 

E tan~^{i'2qk) 



(C.l) 
(C.2) 

(C.3) 
(C.4) 



dS 



\£\<U!/2 

1 



- E^) V2~q 
+ 8£:2 2u''S'' 




tan ^{v'lQ.k) 



Cj\£^ - £i) \ 2E E^ 
Here /< is real, but /> has both real and imaginary parts. We now calculate 



(C.5) 



Re[l>] = ^ / d£ 



+ 



F Jec>\S\>iij/2 

1 

f dS 



E tanh ^{I'sqk) 



2c2- 



^ep J\e\>e. 



2E 
E 



E^ 



k_ 



tanh ^{v^qk) 



tanh~^[l/(i/35^)] 



^sq 



+ 8S- 



k taiilr'fl/ (//,/]/,•)] 



where Sc is the energy £ which is the root of the equation u^qk = 1. In the weak-coupling 
limit, we can carry out the energy integration over the Fermi surface. In that case, the 
integration over the energy \£\ > is negligible and £c may be assumed to tend to 00. 
Under this condition, putting k c:^ k^ c::^ kp, in the zero temperature limit, carrying out 
the Fermi surface integral, expanding the terms like tan~^(gA;) and tanh~^(gA;) up to the 
second order in q, we obtain 



Re[D+] ~ 



Aep 



dS 



E 



E 



{Qkf, 



2E{S'-Si)^ 
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3 



+ - 



4^^ 



2^2 



(C.7) 



4e, 



2 



3uj 



r,2 



[qk^f f 20AVl^^ - 4A^(1 - 2a;^) sm-^(Lj) 



+ 



(gfc^ /4A 



VT— ciJ^ sin ^ uj 



where lo = u;/(2A). For u « 1, we can expand 

1 



sin ^(Vl — uj"^) — sin ""^(l — -ciJ^^ 



TT 



1 



(C.8) 



(C.9) 



2^2 12 

Similarly, expanding sin~^(ci}) and retaining the terms lower order in u, we can 
approximate the square bracketed part of Eq. (IC.Sp as 

{qkY 



12 



2 + l^' 





2{qk, 



9(2A)2(1 - c^2)3/2 + 9(2A)2 

Furthermore, keeping the terms up to second order in uj and neglecting the product 
terms like [qk^Yuj"^, the above expression reduces to 



_2 _ 4 / 

3 Ua. 



Equating this to zero, we find the root u = -^qk^ which implies 
1 



UJ 



V3 



(C.ll) 



(C.12) 



where Pq = Tiq and = hk^/m. 



